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#Include <stdio.h> 
#include <stdlib.h> 
jfinclude <memory.h> 
#include <math.h> 
jfinclude <time.h> 

#define MAX_VALS 20000000 // Maximum number of values to be sorted 

#define MASKWIDTH 8 // Width of the mask to use by Linear Sort 

#define MAXCHILDREN 256 // This should be set to 2 A MASK_WIDTH 

#def ine SEED_INCREMENT 473293813 // Used by the random number generator 
fldefine M0D_VAL 10000 // Values to be sorted range 0 - M0D_VAL-1 

typedef struct val_type 
{ struct val_type *next; 
int value; 

} i 

struct val_type *root, initial__data[MAX_VALS] ; 
unsigned long int valuesjnask, startingjnask; 

int num_vals, initialjrightmost, sortedvals [MAX_VALS] , target, cycles; 
clock_t before, after; 



void prepare jiata (void) 
{ struct val_type *tval; 
int i, seed=SEED_INCREMENT%MOD_VAL; 

values_mask=0; 

starting_mask=0; 

cycles=0; 

initial__rightmost=0 ; 
target=0; 

// set up the values to be sorted 

root=NULL; 

values_mask=0; 

for (i=0; i<num vals; i++) 

{ tval=&(initiaT_data[i] ) ; 

tval->next=root; 

tval->value=seed; 

values_mask=values jnask j seed; 

seed= (seed+SEEDJNCREMENT) IM0D_VAL; 

root=tval; 



for(i=0, start ing_mask=0; i<MASK_WIDTH; i++) // Build the mask 
{ starting_mask=starting_mask*2+l; } 

for (initial rightmost=l; starting_mask<values_mask; ) // find masking start 
{ initial _rTghtmost++ ; 
startingjnask* =2; 



FIG. 7 A 
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void linear_sort (struct val type *curr, int count, unsigned long int 
mask, int shift, Int rightmost) 



{ int i, c, t, children_count[MAX_CHILDREN] ; 
struct val_type *tval, * children [MAX_CHILDREN] ; 

if ((count<=l) | | (mask<=0)) 
{ for (i=0; i<count; i++) 
{ sortedvals [target] =curr->value; 
target++; 



} 



} 

return; • 
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memset(& (children) , 0, sizeof (children) ) ; 

memset (&(children_count) , 0, sizeof (children_count) ) ; 



for (c=0; c<count; C++) 

{ i=(curr->value & mask) >> (rightmost- 1) ; 

tval=curr; 

curr=tval->next; 

tval->next=children[i] ; 

children [i]=tval; 

children_count [i] ++; 



mask=mask>>shift; 
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rightmost-=shift; 



for (c=0; C<MAX CHILDREN; C++) 
{ if (children [c] ) 

{ linear_sort (children [c] , children_count [c] , mask, shift, rightmost); } 



FIG. 7B 
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void quicksort (int loO, int hiO) 
{ int lo = loO; 

int hi = hiO; 

int pivot, t; 

if (lo >= hi) { return; } 
else if( lo == hi - 1 ) 
{ if (sortedvals [lo] > sortedvals [hi] ) 
{ t = sortedvals [lo] ; 

sortedvals [lo] = sortedvals [hi] ; 

sortedvals [hi] = t; 

return; 



pivot = sortedvals [ (lo + hi) / 2] ; 
sortedvals [ (lo + hi) / 2] = sortedvals [hi] ; 
sortedvals [hi] = pivot; 

while ( lo < hi ) 

{ while ( (sortedvals [lo] <= pivot) && (lo < hi)) 
{ lo++; } 

while ((pivot <= sortedvals [hi] ) && (lo < hi ) ) 
{ hi--; T 

if (lo < hi) 

{ t = sortedvals [lo] ; 

sortedvals [lo] = sortedvals [hi] ; 

sortedvals [hi] = t; 



sortedvals [hiO] = sortedvals [hi] ; 
sortedvals [hi] = pivot; 
quicksort ( loO , lo- 1 ) ; 
quicksort (hi+1, hiO) ; 



FIG. 7C 
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void main (void) 

printf ( "#_Values\t\tLinear\t\t\tQuicksort\n n ) ; 

for (num_vals=1000000; num_vals<=MAX_VALS; num_vals+=1000000) 
{ prepare_data ( ) ; 
before=clock() ; 

linear_sort (root, num_vals, start ingjnask, MASK_WIDTH, initial_rightmost) ; 
after=clock() ; 

printf ("%10d\t%10d\t%10d\t\ num_vals, cycles, after-before) ; 

build_dataset () ; 
before=clock() ; 
quicksort (0, num_vals) ; 
af ter=clock() ; 

printf ( n %101\t%10d", cycles, after-before); 
printf (»\n») ; 



void build dataset (void) 

{ int i, high, low, avg, counts [MOD_VAL] ; 

cycles =0 ; 

sortedvals [0] =SEED_INCREMENT%MOD_VAL ; 
for (i=l; i<num_vals; i++) 

{ sortedvals [i] = (sortedvals [i-1] +SEED_INCREMENT) %MOD_VAL; } 
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The following source code sample contains both the Linear Sort and the Quicksort Algorithms, 
tfinclude <stdio.h> 
jjinclude <stdlib.h> 
if include <string.h> 
§ include <memory.h> 
ft include <time.n> 
#define MAX VALS 1000000 
ffdefine TEST INCREMENT 10000 
ffdefine MAX STR LEN 20 
ffdefine MAX"CHILDREN 256 



// Maximum number of values to be sorted 
// Maximum number of values to be sorted 
// Maximum length of strings to be sorted 
// 256 because the Mask Width here is 8bits 



typedef struct val type 
{ struct val_type T next; 
char *value; 

}■ 

struct val_type *root, initial_data[MAX_VALS] ; 

long num vals, target, cycles, moves, compares; 

char *so?tedvals[MAX_VALS] , raw_data [MAX_VALS] [MAX_STR_LEN] ; 

clock t before, after; 

FILE Unfile; 

void prepare_data(void) 
{ struct val_type *tval; 
int " i; 

target=0; 

// set up the values to be sorted 
root=NULL; 

for (i=0; i<num vals; i++) 
{ tval=& (initial datafi] ) ; 

tval->next=roof ; 

tval->value=&(raw_data[i] [0]); 

root=tval; 



FIG. 8A 
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void linear sort (struct val type *curr, int count, int level) 
{ int i, c,~t, children_count LMAX_CHILDREN] ; 
struct val_type *tval, Children [MAX_CHILDREN] ; 

if (count==l) 

{ sortedvals [target] =curr- >value ; 
target++; 
return; 

} 

memset (& (children) , 0, sizeof (children) ) ; 

memset (&(children__count) , 0, sizeof (children_count) ) ; 



for (c=0; c<count; C++) 
{ i=curr->value [level] ; 

cycles++; 

if (i==0) 

{ sortedvals [target] =curr->value; 
^ target++; 

else 

{ tval=curr; 
curr=curr->next; 
tval->next=children[i] ; 
children [i]=tval; 
children_count [i] ++; 



} 

for (c=l; c < MAX_CH I LDREN ; C++) 

{ if (children [c] ) 0 
{ linear_sort (children [c] , children_count [c] , level+1); } 

} 



void validate_sort (void) 
{ int i; 

for (i=l; i<num_vals; i++) 

{ if (strcmp (sortedvals [i-1] , sortedvals [i] ) >0) 

{ printf( n sort error=> %d: [%s] [%s] \n" , i, sortedvals [i-1] , sortedvals [i] ) ; 
return; 



printf(" OK 11 ); 



FIG. 8B 
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void quicksort (int loO, int hiO) 
{ int lo = loO; 

int hi = hiO; 

char *pivot, *t; 

if (lo >= hi) { return; } 
else if ( lo == hi - 1 ) 

{ if (strcmp (sortedvals [lo] , sortedvals [hi] ) >0) 
{ t = sortedvals [lo] ; 

sortedvals [lo] = sortedvals [hi] ; 
sortedvals [hi] = t; 

compares++; 
return; 



pivot = sortedvals [ (lo + hi) / 2] ; 
sortedvals [ (lo + hi) / 2] = sortedvals [hi] ; 
sortedvals [hi] = pivot; 

while ( lo < hi ) 

{ while ( (strcmp (sortedvals [lo] , pivot ) <=0) && (lo < hi)) 
{ lo++; 

compares++; 

compare s++; 

while ( (strcmp (pivot, sortedvals [hi] ) < = 0) && (lo < hi ) ) 
{hi--; 

compares++; 
if (lo < hi) 

{ t = sortedvals [lo] ; 

sortedvals [lo] = sortedvals [hi] ; 
sortedvals [hi] = t; 
moves ++ ; 

)' 

sortedvals [hiO] = sortedvals [hi] ; 
sortedvals [hi] = pivot; 
quicksort (loO , lo-l) ; 
quicksort (hi+1, hiO) ; 



FIG. 8C 
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void build_dataset (void) 
{ int i, c=0, m=0, p=0; 

inf ile=f open ("strings. dat", "r") ; 
for (i=0; i<MAX_VALS; i++) 
{ fscanf (infile, n %s\n\ &(raw_data[i] ) ) ; 
if (strlen(raw_data [i] ) >ra) 
{ m=strlen (raw_data [i] ) ; 



fclose (inf ile) ; 
printf ("max string length=%d at %d\n", m, p) ; 



void reset_dataset (void) 
{ int i; 

for (i=0; i<num_vals; i++) 

{ sortedvals [i] =&(raw_data [i] [0] ) ; } 



void dump jiataset (void) 
{ int i; 

for (i=0; i<MAX_VALS; i++) 
f printf (»%d: %s\n", i, rawjiata [i] ) ; } 



{ printf ("%d: %s\n\ i, sortedvals [i] ) ; } 



void main (void) 

{ 

build dataset () ; 

printf ( " \t\tQuicksort\t\t\tLinear\n" ) ; 

printf ("#_Values compares moves clock cycles clock\n n ); 

for (num vals=TEST_INCREMENT ; num_vals<=MAX_VALS; num_vals+=TEST_INCREMENT) 
{ reset__dataset () ; 

compares=0; 

moves=0; 

printf ("%10d ", num_vals) ; 
before=clock() ; 
quicksorts, num_vals-l); 
after=clock() ; 

printf ("%10d %10d %6d", compares, moves, after-before); 

cycles=0; 
prepare_data ( ) ; 
reset_dataset () ; 
before=clock() ; 

linear_sort (root, numjvals, 0); 
after=clock() ; 

printf (" %10d %6d", cycles, after-before); 
printf ("\n") ; 




p=i; 




FIG. 8D 
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